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Abstract 

We investigate self-moving particles which prefer to hop with a certain turning angle 
equally distributed to the right or left. We assume this turning angle distribution to 
be given by a double Gaussian distribution. Based on the model of Active Brownian 
particles and we calculate the diffusion coefficient in dependence on the mean and 
the dispersion of the turning angles. It is shown that bounded distribution of food 
in patches will be optimally consumed by the objects if they hop preferably with a 
given angle and not straight forwardly. 
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1 Introduction 



In 1827, when the British botanist Robert Brown discovered the erratic motion 
of small particles immersed in a liquid, he considered them first as living 
entities. Later on Einstein and Smoluchowski [1] have shown that the behavior 
of Brownian particles are due to the molecular agitation, i.e. collisions from 
the side of the solvent molecules only. 

The most fascinating formulation was given by Paul Langevin [2] who added 
random forces to the equation of motion, in addition to the viscous forces 
acting on the particles in the fluid. Recently, this latter formulation was used 
to study active particles which are self-moving particles thanks to some energy 
supply and, thus, remember the entities as they were considered by Brown. 
This concept of active Brownian particles is easy, negative friction at smaller 
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velocities resembles a pump of energy. This energy is transformed into kinetic 
one up to some velocity (for higher ones, dissipation wins again) and the 
particles starts to be moving along their given direction. In case of additionally 
acting forces, in particular random ones, the individuals change the direction 
as well as the modulus of their velocity. 



Several interesting phenomena were studied in the past. Simple models of ac- 
tive particles were studied already in many earlier works [3-6]. For broader 
reviews one might refer to [7-9]. Other previous versions of Active Brownian 
particle models [10,11] consider more specific activities, such as environmental 
changes and signal-response behavior. In these models, the Active Brownian 
particles (or active walkers, within a discrete approximation) are able to gen- 
erate a self-consistent field, which in turn influences their further movement 
and physical or chemical behavior. This non-linear feedback between the par- 
ticles and the field generated by themselves results in an interactive structure 
formation process on the macroscopic level. Hence, these models have been 
used to simulate a broad variety of pattern formations in complex systems, 
ranging from physical to biological and social systems [10,12,13]. 



The concept of Active Brownian particles was also applied to the motion of 
small water flees, so called Daphnia [9, 14] where it can also explain some 
collective swarming and curling of the animals and in light shafts [15, 16]. 
These Daphnia exhibit another property during their motion as was recently 
detected in measurements [15]. During a small time step they hop with some 
preferred angle to the left or right of their previous motion. Hence the direction 
of motion is due to some given distribution of turning angles which determines 
the preferred direction after the time step. 



In [17] we applied a random walk model to find the effective diffusion coefficient 
of the random motion with preferred turning angles. For this purpose we have 
used and expression derived by Kareiva and Shigesada (see [18] for details) and 
applied it to the experimental observations [15]. Here we develop a continuous 
model which is based on the stochastic equations of motion in two dimensions. 
Again we find a decrease of diffusion with the mean turning angle and a non- 
monotonous dependence on the dispersion of angles. 



In addition we let the particle consume some non-moving food with constant 
rates. If this food is distributed in bounded regions, it would be suggestive 
to hop with a preferred turning angle. That's why the food consumption in a 
given time is maximized as will be shown in the last section. 
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2 Equations of motions of Active Brownian particles 



The heart of Active Brownian particles is an energy pump modeled by negative 
dissipation at small velocities. Thus we suppose that the relaxation coefficient 
of the velocities is negative for small velocities. For large velocities dissipation 
again overwhelms. A good prototypical choice goes back to Rayleigh who 
considered first nonlinear friction as [19] 

7(f) = —a + v 2 (1) 

where v is the modulus of the velocity-vector. This type of friction function 
describes a source of energy if v 2 < a and lowers the velocity in the opposite 
case. This is done upto a stationary velocity which has the absolute value of 
/o3. 



Since Daphnia hop in first approximation in a two dimensional plane we as- 
sume that the position and the velocities are determined by the equations of 
motion: 

dr _ dv 



dt =v, — = -y(v)v + £x0+y/2D v Z(t) (2) 

We assume that r = {ri,r 2 } and v = {i>i,i>2} are two-dimensional vectors, 
respectively, for the position and the velocity of the considered particle. D v 
scales the randomness of the motion which is due to the stochastic sources in 
(2). Therein is a Gaussian stochastic force with unit strength, independent 
components and a 5-correlated time dependence 

(£(*)) = 0, (fc(0fc(O> = (3) 



The new term in (2) is the Lorentz-like force which models the circular motion. 
We assume that 00 is directed perpendicularly to the plane of motion, i.e. 
uj = {0, 0, 00} and a Larmor-like rotation is generated. The sign of 00 determines 
clockwise or counter clockwise motion. Later we will assume that half of a 
population is running clockwise and the second half counter clockwise. 

The situation can be better understood if we change to polar coordinates in 
the velocity space by introducing 

Vl = v (t) cos (0(t)) , v 2 = v(t) sin (</>(t)) (4) 

This change of variables is quite similar to an often used transformation in 
the phase space of nonlinear oscillators [21]. In contrast we consider a trans- 
formation within the two-dimensional velocity space. 

The equations of motion transform into 
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^jj = -l(v) v + cos(» & + sin(0) £ 2 

= w + ^ (- sin(0) £1 + cos(0) £ 2 ) • 

Without noise the velocity vector relaxes to a constant modulus v 
moves with constant angular frequency uj. 



In Stratonovich-interpretation the corresponding Fokker-Planck equation for 
the density in the v, 0-space then follows [20, 22, 23] 



(5a) 
(5b) 
a which 



d_ 

dt 



P(v,<j>,t) 



d_ 

dv 



■y(v) v — 



dP 



d(j) ^ v dv 2 



d 2 P D„ d 2 P 



+ 



v 2 d(f) 2 



Its stationary distribution can be easily found 
P (v,(f>) = Nv exp (^-jy Jl{v)vdv^j, 



(6) 



(7) 



which is uniformly distributed in 0. We see that the rotational motion has 
no influence on the stationary density. We mention that the Larmor rotations 
are inducing a circulating probability flow. The most probable values can be 
easily found from (7) 



7(6) - 



D 



0. 



(8) 



In the limit of weak noise and strong pumping D v /a — > we approximate 
the velocity distribution by a 5-function around the most probable values 

(v > 0,v > 0) 

Po(v,<t>) = ^-6(v-v) . (9) 



3 Mean square displacement 



We will characterize the motion of the active particles by its mean square 
displacement. In our case without external confinements the random motion 
induces linear diffusion in the long time limit at large length scales. Hence the 
mean square displacement is determined by the diffusion coefficient which is 
in two dimensions 

((f(t)-f(0)f) = 4D r t. (10) 

For a more detailed discussion of linear diffusion and the connected Einstein 
relations see [24, 25] . 
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In order to calculate D r we integrate the velocity, which yields 

f(t) = f(0)+ f'dhvih). (11) 
J o 

Therein we agreed that the particles have been started at t — in r(0). By 
this way the mean square displacement reads 

'fit) - f(0)) 2 ) = J* dh jf dt 2 {vitjvfo)) , (12) 

and is traced by to the velocity correlation function. Crossing to radial and 
angle variables in velocity space (Eq.(4)) gives 

K(t) = (vihWh)) = (u(fi)u(t 2 ) cos(0(*i) - 0(t 2 )) . (13) 

Our main approximation is a widely used decoupling of v(t) from the dynamics 
of the angle <p{t) [26]. Then for the correlation function we find 

K(t) ~ (u(*iM* 2 )> (cos(0(*O - 0(t 2 ))> (14) 

More restrictively we require that the modulus of the velocities is approxi- 
mately fixed which results in (v (ti)v (t 2 )) = v 2 and v is due to Eq. (8). 

Hence the mean square displacement is due to the correlations of angles of 
the velocity vector. From Eq. (6) we get for the angle-distribution the Fokker- 
Planck equation 

d dP d 2 f 



with the diffusion coefficient = D v /v 2 . Solving (15) with t 1 > t 2 gives 

(01 - u(ti - t 2 ) ~ <f> 2 ) 



P(<f>i,h\<h,h) = i ^=exp 



(16) 



This expression is used for averaging the cosine of the increment in angles 
<f>(h) — <P(h)- Afterwards the time integrals can be calculated by standard 
calculations and a lengthy expression is obtained. We give here the long time 
asymptotic behavior where a linear increase becomes dominant over decaying 
and constant items. It has a simple shape as 

f(t) -f(0)) 2 ) = 2 v 2 D f**^ 2 t + O (const, cos(cjt), exp(-D^)) (17) 
Hence the diffusion coefficient becomes 

Dr = {tt^wd 2 ) • (18) 
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The first factor agrees with an expression derived by M ikhailov and Meinkohn 
[27] in case there is no preference in the turning angle {u = 0). The second 
factor in the brackets decreases in case of uo ^ 0. This solution for the spatial 
diffusion coefficient is in qualitative agreement with previous results based on 
a random walk model [17]. 

We point out that the diffusion coefficient does not depend on the sign of u, 
hence not on the direction of rotations. Therefore the derived expression is 
true also if part of the population hops clockwise and the other part counter- 
clockwise. 

The general dependence of the diffusion coefficient (18) on noise D v and mean 
turning angle is shown in Figure 1. It decreases monotonically with the mean 
turning angle u. Remarkable one finds a maximum with respect to the noise. 
It is located at noise values 

and the spatial diffusion coefficient equals v 2 /4u. 




(a) (b) 

Fig. 1. Spatial Diffusion D r depending on external noise D v and the preferred turn- 
ing angle to (a) or the velocity v in the stationary case (b) 



4 Optimal food consumption of finite food patches 



In this section we want to develop a simple approach how much food a single 
Daphnia could consume while it is on its foraging path. Assuming that the den- 
sity of the animal g{f , t) obeys a diffusion equation with the afore-calculated 
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diffusion coefficient D r 



In two dimensions the resulting probability density in space is 
I 



4tt D r t 



cxp 



4£U 



(20) 



(21) 



We assume that the particle during its random motion consumes with constant 
rate k food C which is given by its density c(r , t) . If the food is un-movable 
in space, the consumption of food is described by the dynamics 



d 

— c(f, t) = -k c(r, t) g(f, t) 



(22) 



The latter equation can be solved exactly if the solution (21) is inserted. Simple 
quadrature gives 



c(r,t) = c (f) exp 



- /' kg(f,t')dt' 

Jt 



(23) 



With the definition of the exponential integral [28], Ei(a) = J a °° ^j-dt, we 
obtain in compact shape 



c(r, t) = c (r) exp 



where we have put t = 0. 



k 



4ttA 



■Ex 



(24) 



Let us look at bounded distributed food patches now. As a prototypical dis- 
tribution we take spheres with radius R where C is present. Outside of the 
circle no food is located. For simplicity we assume the circle positioned at the 
origin, i.e. 



co(r) = c if \r\ < R. 



(25) 



and vanishing, elsewhere. Then after a fixed time T we determine the overall 
food which is still left. It results from integration 



C(T) = 2ttc / exp 
Jo 



4ttA 



Ei 



AD r T 



r dr. 



(26) 



Inspection of this expression shows that there exists a minimum with respect 
to the diffusion coefficient D r . Since D r depends monotonously on the mean 
turning angle u of the zig-zag hopping it implies a minimum with respect to 
the mean angle, too. Thus, the variation of the mean angle with change of the 
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effective diffusion coefficient gives rise to a change of the food consumption. 
One might speculate that nature has found by this way the reason to optimize 
the zig-zag hopping. 

The consumed food as a function of the spatial diffusion coefficient D r is shown 
in Fig. 2. A simple explanation can be given as follows. If particles hop straight 
forwardly they will leave the food region, going forwardly and backwardly is 
also of disadvantage. Hence there is a optimal angle which the particle spends 
most time in the food patch with. It does not matter whether they start at 
the center of the food or not. 




5e-06 le-05 1.5e-05 2e-05 

Fig. 2. Consumed food depending on the diffusion coefficient. 



5 Conclusions 

This paper extends the theory of Active Brownian particles with terms which 
obey the motion of individuals with preferred turning angle. This motion could 
have been observed, for example, in zooplankton populations like Daphnia [15]. 
The preferred turning angle distribution, which could have been easily put into 
a Random Walk model [17], comes into play as a Lorentz-like force within the 
Langevin equation. This translation of the preferred turning angle distribution 
into the concept of stochastic differential equations makes the treatment of 
such biological systems more handy in terms of physical interpretations of 
the influencing forces. With the help of the, now more generalized, theory 
of Active Brownian particles a broader variety of swarming systems can be 
described via physical models. 

We calculated the spatial diffusion coefficient, and what is astonishing, is a 
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maximum of spatial diffusion with the noise strength for preferred turning 
angles — n/2 < u < tt/2. The diffusion law (see Eq. (18)) is an extension of 
the theory presented in [27]. This result can be concluded as follows: With 
a preferred turning angle a individual can cover a larger area within a given 
time. This would be an evolutionary advantage as could have been seen in 
Sec. 4. Here we develop a theory of food consumption within a bounded food 
patch. Assuming a Gaussian distributed density for a single individual, and 
using the results from the sections before, a minimum in the left over food 
with the spatial diffusion coefficient has been found. Hence, to cover the max- 
imum of a given area going zig-zag is many times more beneficial than just 
normal diffusion. Within a given time a individual can cover a bigger area 
with preferred turning angle and therefore consume more food, which is of 
evolutionary advantage. 
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